Electronic states and self-doping at a 45° YBa2Cu307 grain boundary 
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The charge redistribution at grain boundaries determines the applicabiUty of high-Tc supercon- 
ductors in electronic devices, because the transport across the grains can be hindered considerably. 
. . , We investigate the local charge transfer and the modification of the electronic states in the vicinity 

0\ • of the grain-grain interface by first principles calculations for a (normal-state) 45°-tilted [001] grain 

I boundary in YBa2Cu307. Our results explain the suppressed interface transport and the influence 

, of grain boundary doping in a quantitative manner, in accordance with the experimental situation. 

The charge redistribution is found to be strongly inhomogeneous, which has a substantial efi^ect on 
transport properties since it gives rise to a self-doping of 0.10 ± 0.02 holes per Cu atom. 

PACS numbers: 73.20.-r,74.25. Jb,85.25.Am 
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Electronic transport in wires and tapes from high-Tc 
compounds is seriously suppressed by structural defects 
and interfaces. In particular, grain boundaries have been 
identified as the main limiting factor, which determines 
the critical current in bulk samples [l| . Grain boundaries 
in YBa2Cu307_5 (YBCO) seem to be depleted of charge 
carriers 0, Q • For this reason, the supercurrent density 
can be enhanced by locally overdoping the superconduc- 
tor (by Ca substitution) close to grain boundaries while 
keeping the grains optimally doped ^4j]. The overdoping 
prevents a reduction of the critical temperature [1] , which 
indicates that suppression of the critical current at grain 
boundaries is connected to interface charging and/or dis- 
tortions of the local electronic band structure. The latter 
would also explain that a perovskite superconductor be- 
haves different from conventional or intermetallic super- 
conductors, like MgB2, where the supercurrent density is 
hardly affected Yet, there is evidence for Josephson 
junction type weak link grain boundaries in MgB2 [7]. 

In perovskite superconductors, the critical current can 
drop down to 1% of the bulk value for a misorientation 
ang le of 10° [1, i For technological applications, 

optimization of low angle grain boundaries is manda tory 
and, therefore, addressed by various groups [HI, [l3, [Ill- 
Large angle grain boundaries, on the other hand, are 
of practical interest as they form excellent Josephson 
junctions. Using bicrystal substrates or the biepitaxial 
principle, based on two different substrates, 45° tilted 
001] grain boundaries can be grown with a high accuracy 
[ill . Because substrate effects can be eliminated from 
transport measurements [l^, such Josephson junctions 
are widely used to gain insight into the superconduct- 
ing state (e.g., the properties and symmetry of the gap 
function, s-wave vs. d-wave) [HI, Il7|. In addition, meso- 
scopic fluctuations of the magnetoconductance point to 
a coexistence of supercurrent and quasiparticle current 

It is surprising that a first principles electronic struc- 
ture calculation, which can account for structural details 
at grain boundaries, is missing so far. Probably, this fact 



is connected to a tremendous demand on CPU time. For 
the simpler Ti-based perovskite insulators SrTiOa and 
BaTiOa, two common substrates for YBCO, analogous 
studies of grain boundaries have been performed by Klie 
et al. [T9I , using the Thomas-Fermi screening approach, 
and by Imaeda et al. j20| . using density functional the- 
ory (DFT). Moreover, superconductor-metal interfaces, 
which are closely related to grain boundaries in high-Tc 
materials [2T','22'], have been addressed by DFT calcula- 
tions. For YBCO-Pd interfaces, as a prototypical exam- 
ple, a net charge transfer of '--^ 0.1 holesper Cu atom off 
the Cu02 planes has been established [23], showing only 
a weak dependance on the interface orientation 24?] . 

Contrary to conventional superconductors, bending of 
the band structure by a variation in the charge distribu- 
tion is strong enough to control the transport in high-Tc 
compounds [ij] because of large dielectric constants and 
small carrier densities (25j . The Thomas-Fermi screening 
length (over which a band-bending is effective) therefore 
grows to the order of magnitude of the superconducting 




FIG. 1: (Color online) Structural setup used for simulating a 
45° YBCO grain boundary (wiggly line), in a projection along 
the YBCO c-ax;is. Large mint-green spheres represent Y/Ba, 
medium blue spheres Cu, and small red spheres O. Whereas 
the square highlights an YBCO unit cell, the outer rectangle 
indicates the boundaries of the supercell. 
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FIG. 2: (Color online) Left: partial Cu 3d and O 2p DOS for the entire supercell, normalized by the number of Cu/0 atoms. 
Center/right: site-projected DOS of selected sites in the 45° YBCO grain boundary: bulk-like coordinated O (center; showing 
two examples) and Cu at the grain boundary compared to bulk-like coordinated Cu (right; both on CuO chain sites). 



coherence length. Consequently, the technical optimiza- 
tion of grain boundaries calls for further insight into the 
local electronic structure. This issue is addressed in the 
present letter for a 45°-tilted [001] grain boundary. Since 
the electronic states are expected to depend strongly on 
the local atomic environment, we use the DFT approach 
to obtain and study a fully relaxed YBCO grain bound- 
ary in a supercell of altogether 317 atoms. Our data pro- 
vide evidence for a characteristic charge transfer between 
the Cu02 chains and CuO chains. We explain this effect 
in terms of the modifications of the chemical bonding. 

Our results are based on density functional theory and 
the generalized gradient approximation, as implemented 
in the WIEN2/c package 26]. By its all-electron scheme, 
this full-potential linearized augmcnted-plane-wave code 
is suitable for dealing with structural relaxation and the 
induced charge redistribution in complex geometries, ap- 
plying a supercell approach [2?'] . From a technical point 
of view, the electronic states at YBCO grain boundaries 
are accessible to the supercell approach, based on three- 
dimensional periodic boundary conditions, since relevant 
effects are proposed to take place on the length scale 
of the lattice constant. Typical experimental screening 
lengths for doped YBCO systems are ~ 4.5 A [H, which 
is an upper limit for the system under investigation. A 
supercell therefore can cover the essentials of the elec- 
tronic interaction, if density functional theory predicts 
the correct screening length. 

The supercell which we use for simulating a 45°-tilted 
[001] grain boundary is displayed in Fig. [l] in a projec- 
tion parallel to the c-axis of the parent YBCO unit cell 
(black square) . Our structural setup reduces the number 
of atomic sites entering the calculation to the necessary 
minimum. A further reduction is not possible since one 
would not maintain both Cu and O atoms in a bulk-like 
coordination. We mention that a (minor) lattice strain is 
present in our supercell by construction. It results from 
a lattice mismatch of less than 1% of the YBCO lattice 
constants between the neighbouring grains in our super- 
cell: five 45°-tilted and seven non-tilted YBCO unit cells 
form the grain boundary, compare Fig.[TJ The strain has 



to be accepted in order to keep periodic boundary condi- 
tions. However, the mismatch is small enough to ensure 
that there is no drawback on the validity of our results. 
For sampling the Brillouin zone a 5 x 5 x 2 mesh with 25 
points in the irreducible wedge is applied. We use 2502 
local orbitals and « 840000 plane waves to represent the 
charge density. The cutoff is set to RKmax — 4. 

Turning to the results of our electronic structure cal- 
culations for the YBCO grain boundary, we give in Fig. 
[5] an overview of the valence states, consisting of the Cu 
"id and O 2p orbitals. The DOS curves which we present 
here and in the following are normalized by the number 
of contributing atoms in order to simplify a comparison. 
The Cu 3c? DOS exhibits a broad structure with one pro- 
nounced maximum centered at about —2 eV. This shape 
refiects the YBCO bulk data except for the fact that the 
width of the —2 eV peak is increased, which is expected 
since the different Cu atoms at a disordered grain bound- 
ary should have slightly different chemical surroundings. 
Likewise similar to bulk YBCO, there are substantial O 
2p contributions in the entire Cu 3d energy range up to 
the Fermi level due to a strong Cu-0 hybridization. The 
latter also explains the Cu 3d admixtures at low energy. 
However, at this point we encounter an important differ- 
ence to bulk YBCO ^] : A distinct portion of electronic 
states appears below —7 eV. These states, of course, are 
related to the modified bonding at the grain boundary. 

In order to analyze the Cu-0 bonding in more detail, 
the DOS can be projected on the atomic sites of interest. 
Corresponding results are given in Fig. [2l In the central 
panel, we compare the O 2p DOS of two (representative) 
bulk-like O atoms, one from a BaO layer and one from a 
CuO chain. The term "bulk-like" is used to indicate that 
an atom is coordinated as in bulk YBCO, i.e. its first co- 
ordination sphere is not affected by the grain boundary. 
For all bulk- like O atoms we obtain a similar behaviour: 
Above some —5 eV the DOS resembles the bulk YBCO 
DOS, whereas two additional peaks appear at lower en- 
ergy. These states trace back to additional Cu-0 bonds 
created at the grain boundary due to the broken transla- 
tional symmetry, see Fig. [T] Yet, because they reappear 
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FIG. 3: (Color online) Partial (site-projected) DOS of selected bulk-like coordinated atoms in the 45° YBCO grain boundary, 
as compared to the corresponding YBCO bulk DOS: Cu in the Cu02 planes (left top), Cu in the CuO chains (left bottom), Y 
(right top), and Ba (right bottom). 



for bulk-like atoms they cannot be described in terms of 
localized states at the grain-grain interface. A possible 
explanation would be a substantial delocalization of the 
hybridized Cu-0 states, which could reach deep into the 
grain. However, our minimal supercell does not allow us 
to resolve the latter. In contrast to the naive expectation 
for d electron orbitals, such a delocalization seems to be 
typical for metallic YBCO 30] . This line of reasoning is 
supported by the Cu 3d DOS. On the right hand side of 
Fig. [2] we show, as a representative example, data for Cu 
atoms from a CuO chain. Both for the bulk-like and the 
boundary atom the new low energy states are present. 

It is well-established that the electronic characteristics 
of copper oxides depend severely on details of the chemi- 
cal bonding and on the doping [3lU32|. Here, our results 
show an interdependence between the modified chemical 
bonding at the grain boundary and the charge distribu- 
tion between the structural building blocks, in particular 
the Cu02 planes and CuO chains. While the creation of 
additional CuO bonds predominantly affects the low en- 
ergy range of the valence states, see the previous discus- 
sion. Fig. [3] demonstrates that there are likewise distinct 
alterations at the Fermi energy {Ep)- For selected bulk- 
like atoms in our supercell, i.e. Cu in the Cu02 planes 
(left top) and in the CuO chains (left bottom), Y (right 
top), and Ba (right bottom), we study the site-projected 
DOS for the grain boundary (dashed line) in comparison 
to the respective bulk YBCO results (solid line). Around 



Ep we find significant shifts of the energetic positions of 
the band edges, indicated by arrows in Fig. [3) 

For the Cu02 planes, the electronic states shift to the 
low energy side. The electron count hence increases near 
the grain boundary, i.e. the hole count is reduced. Quan- 
titative evaluation yields a reduction of 0.10 ± 0.02 holes 
per Cu atom, which is a substantial modification of the 
local doping state and should be accompanied by strong 
changes in the electronic properties of the C uO? planes. 
With reference to the YBCO phase diagram [33], it is to 
be expected that this hole undcrdoping will prohibit the 
transition of the grain boundary into a superconducting 
state. The shift of the electronic states can also be inter- 
preted as a compactification of the Cu 3d DOS for sites 
at a grain boundary, see the left top panel of Fig. [31 This 
corresponds to a reduced delocalization of the hybridized 
Cu-0 states, tracing back to the broken lattice periodic- 
ity. Hence, we attribute the charge transfer towards the 
Cu02 planes to an enhanced spatial localization of their 
nesLT-Ep states. Finally, the additional charge originates 
from the CuO chains, see the left bottom panel of Fig. [3l 
since here the band edges shift in the opposite direction, 
i.e. to higher energy. Due to the small portion of Y and 
Ba contributions to the valence states, these sites play a 
minor role for the redistribution. The hole density in the 
Cu02 planes is shown in Fig. [H by contour plots for the 
area of the square in Fig.[T] The observed hole reduction 
affects almost the entire Cu02 plane to a similar degree. 
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FIG. 4: (Color online) Contour plots of the hole density, i.e. 
the integrated total DOS between the Fermi level and 1.9 eV, 
in the Cu02 planes of bulk YBCO (left) and of the supercell 
(right). The display area is the square shown in Fig. [T] with 
Cu in the corners and O in the centers of the edges. 

where qualitative deviations of the contours appear only 
perpendicular to the grain boundary. 

In conclusion, we have studied a 45°-tilted [001] grain 
boundary in YBCO by density functional theory and per- 
formed a full structure optimization. As to be expected, 
we obtain significant modifications of the Cu-0 chemical 
bonding due to additional combinations of effective over- 
lap between Cu 3d and O 2p orbitals at the grain-grain 
interface. Beyond, the bulk-like Cu-0 states are subject 
to enhanced localization near the grain boundary as the 
crystal lattice is perturbed. This results in an energetic 



shift of the electronic states in the vicinity of the Fermi 
level, which has almost the same amplitude in the Cu02 
planes and the CuO chains, but the opposite sign. 

As a consequence, there is a substantial charge trans- 
fer between these structural building blocks. The Cu02 
planes are affected by a charge carrier (hole) depletion of 
0.10 ±0.02 holes per Cu atom. A grain boundary, hence, 
can be interpreted in terms of (intrinsically) underdoped 
Cu02 planes, which explains a suppressed transition into 
a superconducting state. In addition, it explains the fact 
that the supercurrent density can be enhanced by a local 
hole overdoping of the grain boundaries [1, Q , since this 
compensates the intrinsic underdoping. Assuming that a 
third of the additional holes enters the Cu02 planes, the 
maximum of the enhancement for 30% Ca doping agrees 
even quantitatively with the obtained first-principles re- 
sults. Finally, going beyond experimental data the 
charge redistribution is foimd to be inhomogeneous with 
depletion and accumulation zones at the grain boundary, 
while there is little charge transfer to/off the grain. New 
experiments are indicated to confirm this observation. 
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